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PATENT 
Case No. NQQ80US 



IN THE UNITED STATES PATENT AND TRADEMARK OFFICE 

Appl. No. : 09/729,939 

Applicant : Rajashri Joshi et al. 

Filed : December 5, 2000 

Titled : Method and System for Representation of Geographical Features 

in a Computer-based System 



DECLARATION UNDER 37 CFR 1,131 

The undersigned, OLE HENRY DORUM, hereby declares that: 

1. I am a co-inventor of the invention described and claimed in the above-identified 
patent application. 

2. On July 10, 2006, 1 signed a DECLARATION UNDER 37 CFR 1.131 along with co- 
inventors, RAJASHRI JOSHI and VUAYA ISRANI. In that DECLARATION we indicated that 
before August 25, 2000, we invented a new method for representing geographic features. Part of 
this new method included, fitting a polynomial spline to the a geographic feature by applying a 
least squares approximation to data points specifying latitude and longitude coordinates to 
generate a plurality of control points for the polynomial spline. 

3. Before August 25, 2000, a source code was developed; a redacted copy of this source 
code is attached hereto (Exhibit 1). The source code fits a polynomial spline to a geographic 
feature by applying a least squares approximation to data points specifying latitude and longitude 
coordinates to generate a plurality of control points for the polynomial spline. Specifically, the 
source code includes a function named 'uniformbsplinelsfit' that fits a polynomial spline to data 
points by applying a least square approximation to generate control points for the polynomial 
spline. Additionally, the function 'uniformbspline* plots the generated control points and the 
polynomial spline line. 
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4. Before August 25, 2000, the source code was tested and found to work. Specifically, 
the source code was tested and provided the output plot attached hereto (Exhibit 2). The blue 
points in the output plot are shape points specifying latitude and longitude coordinates of points 
along a road segment. The control points of the polynomial spline computed by the source code 
are shown as black circles in the output plot and the polynomial spline line is shown in red. 

5. Exhibit 3 is a table that maps the claim elements of each independent claim onto the 
lines of source code of Exhibit 1 and features of the output plot of Exhibit 2. 

6. All statements made herein of my own knowledge are true and that all statements 
made on information and belief are believed to be true, and further these statements are made 
with the knowledge that willful false statements and the like so made are punishable by fine or 
imprisonment, or both, under Section 1001 of Title 18 of the United States Code, and that such 
willful statements may jeopardize the validity of the application or any patent issuing thereon. 





OLE HENRY DORUM 



Date 



NAVTEQ North America, LLC 
222 Merchandise Mart Plaza, Suite 900 
Chicago, IL 60654 
(312) 894-7000 x7365 
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EXHIBIT 1 



Source Code 
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function [Px, Py, Kx, Ky] = 

uniformbsplinelsf it (x,y,xd,yd, t,C,N, rf ,ml,mc, w) 

% Draws a Is fit uniform B-spline curve given the data points, 
% tangent vectors at the endpoints, parameter vector, knot vector, 
% number of control points, number of data points and 
% regularization factor and tangent weight 

for i = 1:C 

for m = 1:N 

Bv(i,m) = Bu(i-1, t(m) ) ; 

end; 
end; 



Bvd(l:4,l) = [-0.5 
Bvd(C-3:C,C) = [ 0 
NodeWeight = 10000; 



-0.5 



0.5 0]'; 
0 0.5]'; 



zeros(Cl); Sy = zeros (C,l); 



% Compute S, Sx, Sy 
S = zeros (C,C) ; Sx 
for 1 = 1:C 
for i = 1:C 

for m = 1:N 

S(l,i) = S(l,i) + 

end 

% Add extra weight to nodes % 
S(l,i) = S(l,i) + NodeWeight 
S(l,i) = S(l,i) + NodeWeight ' 

% Set up least squares equation matrix 
S(l,i) = S(l,i) + ... 



end 

end 

for 1 = 1:C 

for m = 1:N 

Sx(l) = Sx(l) 
Sy(l) = Sy(l) 

end 

% Add extra weight to nodes % 
Sx(l) = Sx{l) + NodeWeight 
Sy(l) = Sy<l) + NodeWeight 
Sx(l) = Sx(l) + NodeWeight 

Sy(l) = Sy(l) + NodeWeight 

% Set up least squ ares equation data matrix 
Sx(l) = Sx<l) + w*| 
Sy(l) = Sy(l) + w*| 

end 





xc = inv(S)*Sx; yc = inv(S)*Sy; 
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1 [Px, Py, Kx, Ky] = unif ormbspline (xc , yc ) ; 

2 function [Px, Py, Kx, Ky] = unif ormbspline (x, y) 

3 % Evaluates the uniform cubic B-spline curve given the 

4 % control points 

6 MB - (1/6) * [ -1 3-3 1; 

7 3-6 3 0; 

8 -3 0 3 0; 

9 14 10]; 

10 

11 MHA = [ 2 -2 1 1; 

12 -3 3 -2 -1; 

13 ooio 

14 1 0 0 0]; 

15 

16 MHB = (l/6)*[0 0 -1 1; 

17 0 0 3 0; 

18 -6 6 -2 -1; 

19 6 0 0 0]; 

20 

21 C = length(x); Delta = 0.1; InvDelta = 1/Delta; 
22 

23 Px = [ ] ; Py = [ ] ; Kx - [ ] ; Ky = [ ] ; GA = [ ] ; GB = [ ] ; 

24 

25 for seg=l:C-3 

26 if (seg == C-3) i = 0:Delta:l; 

27 else i = 0:Delta:l-Delta; 

28 end; 

29 xt = x (seg: seg+3) ; yt = y ( seg : seg+3 ) ; 

30 Pxi=polyval (MB*xt, i) ; Pyi=polyval (MB*yt , i) ; 

31 Kx = [Kx; Pxi(l)]; Ky = [Ky; Pyi(l)]; 

32 Px = [Px; Pxi'J; Py = [Py; Pyi']; 

33 

34 GHA = inv(MHA) *MB* [xt yt] ; 

35 GA = [GA; GHA] ; 

36 GHB = inv(MHB) *MB* [xt yt] ; 

37 GB = [GB; GHB] ; 

38 end; 
39 

40 Kx = [Kx; Pxi(end)]; Ky = [Ky; Pyi (end)]; 
41 

42 % Find cubic spline coefficients, type A 

43 if 0 

44 x = GA ( : , 1 ) ; y = GA ( : ,2) ; 

45 C = length(x); Delta = 0.1; InvDelta = 1/Delta; 

46 for seg=l:4:C 

47 xt = x (seg: seg+3 ) ; yt = y (seg: seg+3 ) ; 

48 Pxi-polyval (MHA*xt , i ) ; Pyi=polyval (MHA*yt , i ) ; 

49 plot(Pxi, Pyi, 'g-'); 

50 end; 

5 1 end ; 

52 

53 % Find cubic spline coefficients, type B 

54 if 0 

55 x = GB(:,1) ; y = GB{:,2) ; 

56 C = length(x); Delta = 0.1; InvDelta = 1/Delta; 

57 for seg=l:4:C 

58 xt = x (seg: seg+3 ) ; yt = y (seg: seg+3) ; 

59 Pxi=polyval (MHB*xt , i) ; Pyi-polyval (MHB*yt , i ) ; 

60 plot(Pxi, Pyi, 'g-'); plot(xtd), yt(l), 'go*); 

6 1 end ; 

62 end 
63 
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EXHIBIT 2 



Sample Output from Source Code 
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